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ABSTRACT 

There is a need for a fast method of analyzing bathythermograph 
traces and this need is approached by the means of a high-speed digital 
computer. The theory behind the computer program is outlined. Both 
synthetic and real data cases are run as examples using both data card 
decks and magnetic tape inputs. The program has been designed to read 
digitized bathythermograph traces and then analyze them objectively by 
Gaussian and non-Gaussian methods for the top, center, and base of the 
main thermocline. Additionally, such features as multiple thermoclines, 
inversions, and thermal transients are identified also and their key 


points are included in the informational data printout. 
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is Introduction. 

Since the introduction of the bathythermograph by Rossby and 
Montgomery (1934), and its further modification and perfection by 
Spilhaus (1938), there has been almost continuous thermal structure 
information funneled into the various collection agencies which has 
been handled manually under procedures outlined by LaFond (1951). 

This data collection has resulted in vast numbers of bathythermograph 
slides that have been noted, and then stored away for future analysis. 
In more recent years bathythermograph data has been stored on aperture 
cards, using techniques similar to those suggested by Sauer (1964). 

The Canadian Océanographic Data Center, Scripps Institution of Oceano- 
graphy, and the Fleet Numerical Weather Central, Monterey, California, 
have further modified this procedure by a method similar to that sug- 
gested by Sauer and Hope (1967), where a device has been invented which 
will take incoming bathythermograph slides and convert them directly to 
digitized data format. This system described by Sauer and Hope (1967) 
consists of taking the grid/slide combination and transferring it to an 
aperture card by means of the BT Processor, a modified 3M Processor 
Camera. It has the advantage of preserving the original analog trace. 
The aperture card is then fed through a hopper device directly into a 
Caps Jeffreé projector with a D-MAC curve-follower unit. The time re- 
quired for this unit to process one bathythermogram trace is approxi- 
mately 45 seconds with a claimed accuracy of 1.0% of the full depth 
range and 0.1 degrees Farenheit. After exhaustive testing at the 

Naval Hydrographic Office as described by Stewart (1963) Sauer and 

Hope (1967) of the Canadian Oceanographic Data Center felt this degree 


of reliability in fact was not approachable. Accordingly, they 


converted their system to record temperature to a minimum of every 0.1 de- 
grees Celcius and to the nearest whole meter. 

Papers by Tully (1961, 1964) show that there are certain features 
which must be identified in order to classify a bathythermograph trace. 
Tully states that these few parameters can be used to identify the thermal 
structure in any area of the ocean. Among these features is a potentially 
isothermal layer, which may or may not contain transients. There is a 
seasonal thermocline, in which the gradients may be regular or irregu- 
lar. Below this layer is a sub-thermocline layer in which the gradients 
are less than in the thermocline. Superimposed on these features are 
various types of short term fluctuations which possess varying degrees 
of significance, such as temperature inversions. 

Tully (1964) suggests in this same paper that the thermal structure 
can be defined by the following features: 

(1) Surface Temperature. 


(2) Potential Layer Depth. This figure includes the average 
depth to the top of the thermocline, plus any modifications due to inter- 
nal waves, etc. 

(3) Depth to the Bottom of the Thermocline. 

In many oceanic areas this depth has an associated distinct tempera- 
ture minimum and therefore is very definitive. In areas like the 
tropics it is usually not a significant feature. 

(4) Thermal Magnitude of Thermocline. To researchers this 
may be a very significant feature, because it can signify a boundary 
to restrict fish movement, water mass mixing, etc. 

(5) Other Less Definable Features. There are other parameters 
that may be considered, like the location of temperature minimum in a 


thermal inversion, mixed layer depth, etc. These features cannot be 
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fitted easily into a generalized machine processing format for ready 
identification. They are not common features, and may have special 
variables that depend solely upon the area of the ocean being investi- 
gated, or upon the season of the year. 

As Boston (1966) states, "until the information on bathythermograph 
slides can be conveniently analyzed by machine, much information about 
the thermal structure of the upper layers of the ocean is quite likely 
to remain hidden in an ever growing mountain of smoked and/or gold film 
glass slides, photographs of same, or punched cards."' Two distinct 
problems arise: (1) the digitizing of the bathythermograph data to a 
form ie Fen be processed by high-speed computers and data handling 
systems; and (2) specification of a program that will take this digi- 
tized data and process it for use in oceanographic research; the pro- 
gram needs the generality to be able to handle any shape trace, and flexi- 
bility so that further routines may be added or subtracted at will with- 
out destroying the basic integrity of the program. As we have seen, 
some progress is being made in the handling of the first problem. 

The second problem is the topic of this thesis; the objective 
analysis of digitized bathythermogram traces. The standard features 
as suggested by Tully (1964) are identified and printed out in a concise 
informative way. One or two additional features, like the identifica- 
tion of thermal inversions and transients, are included in an attempt 
to show that further features may be added to this program according to 


the requirements of the individual user. 
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Ze Method of Analysis of Thermal Structure. 


There will be two separate analytical methods developed for handling 
the objective analysis of the thermal structure. In his technical paper, 
Boston (1966) suggested the following two methods, on which this thesis 
is based. One method consists of a Gaussian temperature structure being 
assumed for the thermal parameter discussed. The second method is a 
common finite difference method, utilized to find inflection points in 
the various trace slopes. 

A prerequisite to both methods is that the data from the original 
analog bathythermograph trace be reduced to an ordered set of depths 
read at some predetermined temperature interval. From this ordered set, 
the thermal structure can be defined objectively in terms of certain pri- 
mary and secondary features. Both methods require only the first and 
second finite differences with depth; therefore only data from depths 
corresponding to some constant thermal interval (conveniently 0.5 to 1.0 
degrees) are required to handle the analytical technique. However, it is 
convenient to include temperature data corresponding to these points, as 
in the computer output there are locations where information as to the 
thermal magnitudes are as important as the spacial magnitudes. 

Thermal structure, as defined in this thesis, is the relationship of 
the temperature, T, to depth, Z, as a function of time and position. 
Furthering earlier ideas, Tully (1964), Tully and Giovando (1963), and 
Uda (1963) show in their papers that the thermal structure changes with 
season, latitude, and even oceanic region. However, Tully (1964) also 
suggested that the number of different features for a given application 
is small. Section A considers those features common to most oceanic 
regions within the upper 300-400 meters of the ocean. Section B considers 


those features that are of major concern in such fields as antisubmarine 
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warfare, biological forecasting, and fisheries that are not discussed 
or included in Section A. 

A. Primary Features: 

The primary features as suggested by Tully (1964) and Boston (1966) 
consist of several convenient parameters; (1) surface temperature, (2) 
depth of potential layer depth (upper isothermal), (3) depth of the 
bottom of the thermocline, and (4) magnitude of the thermocline. 

The potential mixed layer is that as defined by H. O. Pub SP-105, 
Vol 5, and is considered to be isothermal, although this "isothermal" 
region accepts small positive or negative gradients up to 0.3 degrees 
Farenheit per hundred feet. The termination of this nearly-isothermal 
layer is usually abrupt enough to clearly mark the region where the 
thermocline or other feature takes place. The depth of the mixed layer 
will be defined as a the top of the thermocline. Likewise, Lae the 
base of the thermocline, can be considered the top of the lower nearly- 
isothermal area. 

Once the top and bottom of the thermocline have been defined, with 

Te. as the temperature at an » and Tppas the temperature at Zbpt ; 


the definition of the magnitude of the thermocline is 


LAT = nie mists (1) 


Bre Secondary Features 
The secondary features are usually considered as time varitions in 
the primary features of Section A, and can occur at any point in the 


vertical water column being considered. 


sigs Transients 


The term transient will be defined here as referring to short term 
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fluctuations of the thermal structure. It is usually the result of wind 
mixing, and vanishes over a period of a day or more. The program pro- 
duced for this thesis will recognize these transients, assuming that 
these thermal transients take the usual form of small but sharp gradients 
that appear as irregularities on the main thermal structure. Transients 
are often found in association with the well known "afternoon effect." 
This program will aid in the verification of the prediction methods out- 
lined by Tabata and Giovando (1962). 
2. Multiple Thermoclines 

Heating and wind mixing conditions can also produce more than one 
stable thermocline. The growth of two stable thermoclines has been well 
documented by Tully and Giovando (1963). The occurrence of these multiple 
thermoclines is sufficiently common and vital to underwater acoustic pre- 
diction systems that they merit identification. 

3. Inversions 

Inversions are a common oceanic feature, an example of which has 
been described by Roden (1964), and must be identified if one is to ac- 
curately predict underwater sound transmission. The thermal inversion 
takes form when the temperature below the main thermocline (below Zot) 
decreases with depth to a minimum, then increases markedly, and finally 
decreases again, but now usually at a more reduced rate than the first 
decrease. These thermal features are of extreme interest in underwater 
communications as the point of minimum temperature in this thermal 
structure forms the axis of a "sound channel" formed by the thermal bend- 
ing of acoustic rays by this inversion. 

The following methods will be used to identify these primary and 


secondary features within the computer program designed for this thesis. 
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Cc. The Gaussian Thermocline 


A gaussian, or normal, distribution (see Figure 1) of temperature, 


T, as a function of depth, Z, is given by 


-1/2)2" 
T(Z) = 1. “ ee (2) 


{21r 


whose corresponding frequency function is 


47 (~1/5 \2? 
e.: » ec (3) 
ia {271r 


If the distribution function is T(Zm) » where © greater than O and 
Oo 


m are constants, then the frequency function becomes 
2 


dT(z) e 20% (4) 


=— 
dz oat 


The function described by (4) has the following characteristics: 





(a) it is symmetric about Z =n, 
(b) it has two symmetric points of inflection at 
Z=mto , and 

(c) it has a maximum rate of change at Z = ™M* o|3 
A change in the numerical value of m causes a displacement of the curve 
in the vertical direction, Z, but does not alter the form of the curve. 
A change in &@ has the effect of altering the scale of both coordinate 
axes. The smaller the value of ©O , the more concentrated the curve is 


about the point Z =m. If the water surface is defined as Z = 0, and the 


depth considered to be infinite, then 


ri2m)= to} e FF dz 
& of2T J. 4 


is the function defining the vertical temperature distribution. The 
center of the thermocline, or the center of the curve of the T (4m) 


distribution, is at 
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Z° | = di 
Normalized Temperature dZ 
7 = SSeS Se = Be a ea ---miBo 
dT 
aZ 





Qo. | | | O5 , | | | 10. 
(after Boston 1966) 
Fig.! Gaussian thermocline and its frequency 
function (first derivative) 


L6 


ge (6) 


The top and bottom of the thermocline are then defined as being located 
at the points of the maximum rate of change of slope of the T (43) curve; 


namely, 


Zt =-meol3 (7) 
Zip = -m -oJ3 (8) 


which are illustrated in Figure l. 

The handling of bathythermograph data leading to the application of 
results (6), (7), and (8) is very simple, since values of m ando@ may be 
computed directly from the (t,Z) distribution curve. Consider the first 
two moments of the (T,Z) curve about the origin, Z=0, whereby the (nth) 
moment of the curve about the origin is meant the mean values of Zo 
defined as Za » where n is a positive integer. Then with the normal 


curve, the first moment is identical with the center of the thermocline: 
Tlo) 


a (9) 
oe = T(0) -T(-00) wt al 


where T(O) is the temperature at the surface Z = 0 and T (=-@) is the 
temperature at the bottom of the lower isothermal layer which, in this 
idealized case, is assumed to extend to the bottom of the ocean. The 


second moment 


T(0) 


=2 : | a 10 
é TO) - T(-00) te oT 2 


¢ ? 2 eo e 
can be used to compute @ , since the variance, © , is given by 


Beem 2 oo 
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In practice, several simplifications are possible. The temperature T (-@) 
is measured at the bottom of the bathythermograph cast. The differentials, 
dT, become finite differences, AT, and the integrais are replaced by sum- 
mations. If the depths of the temperatures are read at one degree intervals, 
for example, then the T's can be replaced by unity. 

One must remember that although integration is a continuous process, 
the less accurate summation over finite intervals is not. The larger the 
interval, the less accurate is the summation. The working value of sigma, 
as computed by the routine, cannot exactly equal the correct value obtained 
from integration over the normal curve. Thus equations (7) and (8) must 
be modified by some weighing factor depending upon the temperature inter- 
val chosen. Allowing this factor to be defined as, sk, equations (7) and 


(8) become 
Ze => sklo (7A) 


aa ~m - skle (8A) 


where a temperature interval of 0.5 degrees gives sk equal to 1.47 and a 
temperature interval equal to 1.0 degrees gives a value for sk equal to 
130. 

In summary, the objective definition of the top, bottom and center 
of a Gaussian thermocline requires only readings of (a) surface temperature, 
(b) bottom temperature, and (c) depths at some constant pre-chosen tempera- 
ture interval. 


The center of the Gaussian thermocline is 


are ae. ere 


Ze 2h = D2 
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Since 





| YY 2 
Z(Tn); is i Te Alive (13) 
Then 


ae = Z (in), a has (14) 


The top and bottom of the Gaussian thermocline readily follow from (7A) 
and (8A). 

One problem remains: the decision as to whether or not a normal 
distribution exists; lacking this distribution, the inaccuracies built 
into this process are not tolerable. This determination can be accomp- 
lished by computation of the third moment. However, in the following 


section it is apparent that this computation is not necessary. 


D. The non-gaussian Thermoc line 


The analysis of the more common type of thermal structure, the non- 
gaussian thermocline, requires the computation of finite differences on 


Z(Tn). The first finite difference will be defined as 


AZ(To) = Zh) aii) (15) 


and the second finite difference as 


A? Z(Tn) = Alena) - AZ (Ta) (16) 


Generally, the first finite difference is always positive. The second 

finite difference will be negative until, in this case, the point of 

inflection of the (T,Z) distribution curve is reached, at which point 
AL(Tn) becomes non-negative. If the first non-negative second 


finite difference is AZ (Ta) , then the center of the thermocline is 


alae.) 
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In order to locate Ze only the upper portion of the (T,Z) distribu- 
tion curve is considered, where a (24 Be synthetic lower portion to 
this curve segment is constructed so that together with the upper real 
segment a normal distribution is formed. From this curve, O is com- 
puted by methods described in the previous section and Z is given by 
formula (7A). (see Figure 2b). 

Similarly, only the lower (7 >7Z<) part of the (T,Z) distribution 
curve is utilized in computing Z,, . An upper portion of the curve 
is added to the segment and a new value of oc is computed. Zp is 


then calculated by using formula (8A). (See Figure 2c). 


EZ Irregular Thermal Features. 


i Transients 

Transients, according to Boston (1966), are usually considered 
to be less than one degree in magnitude and might not be seen by a greater 
thermal interval. Thus, it is necessary to have a small grid spacing, 
say 0.1 degrees, to identify these transients. Spilhaus (1938) relates 
that the bathythermograph is seldom more accurate than 0.5 degrees, 
borne out by the testing of Stewart (1963), due to such problems as 
dull stylus, hysteresis, and calibration. Even if the instrument were 
exact, and hysteresis were not present, the thickness of the trace limits 
reading accuracy to about 0.2 degrees. Hence, although the grid spacing 
chosen for this program is unrealistic at 0.1 degrees, it is retained with 
the possibility of more accurate thermal instruments, or more accurate 
bathythermograph readings, in the future. 

If such accuracy were possible (see Figure 3) transients could be 


identified by simply identifying the sign changes of Ajulais. » Lf the 
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first negative value encountered was AZ(Tr) then the beginning of the 
first transient would be at ZT aes ) » Then the end of the first thermal 
transient and the beginning of the second transient is signified by a 
sign change of AZ(Th) from positive to negative. If this occurs at 
AZCTn) then the end of the first transient occurs at Z(Tn) and 
the beginning of the next occurs at Zi Vue . This process will con- 
tinue until the top of the first major thermocline is encountered; 
generally we are interested in only the transients in the upper isothermal 
region. 

20 Multiple thermoclines 

The analysis of a thermal distribution containing more than one 
major thermocline proceeds along the lines just discussed in the pre- 
vious section on transients. In Boston's article (1966) he states that 
a grid spacing of one degree is generally adequate if there are only two 
thermoclines, and a grid spacing of 0.5 is necessary should there be more 
thermoclines present. 

As in the transient analysis the sign changes of AZ(Th ) are 
examined. If the first positive value is O2(To) then the center of 
the first thermocline would be located at Z Twa =e » The center of 
the second thermocline, Loo » 1S indicated when the sign of AZCTh) 
changes from negative to positive. If this occurs at A2CTn) » then, 
as in the first thermocline, the center is located at Lhe = Zor. 

The centers of all subsequent thermoclines are found by continuing this 
process throughout the length of the trace. 

When the centers of the shallowest and deepest thermoclines have been 
identified, the method outlined in Section D, the non-gaussian thermocline 
will be used to determine the top of the shallowest thermocline and the 


bottom of the deepest thermocline. Then considering the fact that Q2(Tn ) 
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is a measure of the slope of the (T,Z) curve, the maxima in AZ (Tr) 
corresponds to the most nearly isothermal conditions. Thus these maxima 
can be used to separate the thermoclines. Once this seperation has been ac- 
complished, the procedure from Section D can be used again to locate the 
top, bottom, and center of each individual thermocline. Considering the 
fact that the thermal structure is assumed to be continuous, with only 
small spacing between consecutive thermoclines, there is little need in 
making the identification any more accurate, or finer, than the subdivi- 
sions that are indicated by the maxima. Figure 4 shows an example of 


real traces in the ocean with this multiple thermocline feature. 


a Inversions 

Large inversions (+ 0.5 degrees Celcius) are usually found 
where water sources of different origin are close together or where 
strong velocity shears are encountered. Tully and Giovando (1963) docu- 
ment this feature. Uda (1963) and Tully and Dodimead (1963) also show 
that in regions such as the subarctic, where there is excess precipita- 
tion, the summer profile may be characterized by a nearly permanent sum- 
mer inversion located at about 100 meters depth. 

The main feature (see Figure 5) of an inversion is the sign 
change of the slope of the (T,Z) distribution curve. The major points of 
interest to the investigator are usually considered to be: 

(a) upper and lower boundries denoted by Zc, and Zc? 
respectively, 

(b) magnitude of the inversion, ZUPP and TUPP, being the 
spacial and thermal magnitudes, respectively, 

(c) minimum temperature, TMIN, 

(d) depth of minimum temperature (depth of axis of sound 


channel), ZMIN. 
24 
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Referring to Figure 5, it should be noted that this system of analysis will 
require a more elaborate notation system, since the curve is no longer 
single valued (that is, more than one depth, Les will have a specified 
temperature). In order to distinguish between these depths the notation 
a: (Tn) will be introduced where the subscript j refers to the order 
in which a multiple occurrence happens from the surface down. Thus, if Tn 
occurs at three depths, the shallowest depth would be indicated as Z,(Tn) , 
the intermediate depth Z>(Tn) » and the deepest as Z3(Tn) . As the 
trace is read vertically from the surface down, the trace appears to be 
single-valued (one depth has one temperature), then double-valued (one 
temperature at three depths), before returning to a single-valued state 
again. Since the finite differences must be computed, its first and 
second finite differences will be defined as OZ; (Tn) and afb: (Tr), 
respectively. 

The upper and lower boundaries of the inversion occur at the first 
double-valued temperature. If this temperature is Tn, then Z¢,= Z,CT ) 
and Lor = Z>(™m) » The last double-valued temperature, or a lower- 
valued single temperature included between the last double-valued tempera- 
tures, is the minimum temperature, TMIN, and the depth of this minimum 
temperature, ZMIN, is defined from Ly (Trnsaage) » The magnitude of the 
thermal inversion, TUPP, is given by the absolute value of the difference 
between the first and last double-valued temperatures. Once these in- 
version features have been identified, pertinent features above and be- 
low the inversion still must be identified by the methods outlined in 


Section D. 
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ss Example of Input 


Two methods of input to this program are data deck and tape deck. 
These two examples illustrate these methods. 
A. Data Deck Input 
In using a data deck input, where an analyzer does not have 
access to the digitizing machines and must key-punch his data deck, it 
is possible to input both positional and curve data. The data deck must 
be made up as follows (see Figure 6); 

(1) The data must be printed from the bathythermograph 
curve, at some pre-chosen thermal interval, i.e., - 0.5 degrees. The 
data are punched on the data aperture cards by punching the temperature, 
followed by its depth as picked off at the chosen delta-T. Each data 
bit is allowed ten spaces on the aperture card and no other data can be 
included in this ten space unit, so that a maximum of four data sets 
can be punched per data card,: (since there are only 80 spaces allowed per 
card). The thermal data is continued through as many aperture cards as 
necessary to insure complete recording of the entire trace. 

(2) Immediately following the completion of the first 
trace, the operator must insert "Control Card 1"', which is a dummy con- 
trol card consisting of "0.0" substituted for the temperature, T, and 
any "non-zero" value substituted for the depth, Z. This dummy card 
signals the main program that the reading routine has completed reading 
the data trace and is ready to read the positional data, if it exists, 
or proceed on with the analysis of the trace. 

(3) Following "Control Card 1" we insert the positional 
data card, if we desire to read in this information. This card contains 


the geographical position, both latitude and longitude to the nearest 
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Trace Data Card 





Control Card 1 
0.0 De oO) 


Positional Data Card 


36° 459.55 Pils 3. au 1 1 23 10 68 231630 639957 





Control Card 2 
O20 0.0 


Control card 2 


Control card 1 (nth) 
Position Card (nth) 
Trace Data cards (nth) 


Position card (lst) 
Control Card 1 (lst) 


Trace data cards (lst) 





Fig.6 Card Format and Data Deck Assembly 
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hundredth of a minute, the local date, the Greenwhich Mean Time and date 
of the bathythermograph cast, and the serial number of the bathythermo- 
graph. If the positional data card is used, then every position must be 
filled in exactly as noted in Figure 6 and Appendix II. If the neces- 
sary data is not known, then the unknown spaces must be filled with "0'''s. 
The format on this card is absolutely vital and inflexible. 

(4) Next we repeat the above procedures (1) through (3) 
for each bathythermograph trace desired. It is not necessary to insert 
a positional card with each trace, but if the positional card is used it 
must be complete, as mentioned above. 

(5) When the entire data package has been made up of trace 
data and positional data cards, then "Control Card 2'' is inserted as the 
last card to signal the program that the data is finished and the pro- 
gram can terminate. "Control Card 2'' consists of a normal data card for- 
mat, except the "0.0'' is inserted in position for both temperature, T, 
and depth, Z. This is always the last card in the data deck (if we do 
not consider the programming control cards used by the system), and will 
always follow the positional data card for the last trace or if this data 
is not used, then it will follow the last "Control Card 1.'' The entire 
data deck is assembled as shown in Figure 6, and may be of any length. 
The only programming prerequisites are that the thermal inverval be a 
constant for each individual trace, although it is not necessary for all 
the traces to have the same interval, and the input format, Figure 6, 
must be followed rigorously. All other necessary computations, such as 
sigma, weighing factor, etc., are contained or computed within the 


program. 
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B. Magnetic Tape Data 


When the input data can be digitized, by such facilities as the Fleet 
Numerical Weather Central, Monterey, or the Canadian Oceanographic Data 
Center, a few modifications must be made to the basic program reading meth- 
ods as they presently exist in Appendix II. In order to use magnetic tape 
input it is necessary to remove the cards that are serialized (right hand 
margin) 0185 - 0204 and insert in their place new cards designed to read 
the magnetic tape. It is important to note here that the format of the 
tapes and even the type of tapes makes a vast difference in the make up 
of the reading cards and great care must be used in inserting the cor- 
rect control and reading cards in order to insure not shutting down the 
computer. Before using tapes, be sure to check the individual instruc- 
tion manual and the individual tape specifications and rigorously follow 
the instructions on the format for the make up of these new cards to be 
inserted. The following information must be considered in each case: 

(a) whether seven or nine track tape, 

(b) data density, 

(c) whether even or odd parity, 

(d) whether the data are translated or not, 

(e) if translated, what machine language, 

(f) record length, 

(g) whether or not the data are blocked, 

(h) if blocked, what size. 
All this information is vital in making up the correct system control 
cards for the tape reading unit. Once the system control cards are cor- 
rectly made up and inserted around the program, then a normal "read" card, 


for example; 
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READ (8,1, END=60) 
1 FORMAT(8A5) 


is used, where, in this example, we are reading on tape reader "8", 


sixty traces consisting of eight data sets per trace. 
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an Computer Program Operation. 


The computer first reads in the data for an individual trace, as 
outlined in Section 3, and then reads in the positional data (if existent 
in data deck). Once the individual trace is read into the program, the 
computer takes over, using the analysis method programmed by this routine 
to start the investigations for the features outlined in Section 2. Be- 
fore starting any further investigations the program looks at the thermal 
interval and choses one of three possible values for the weighing factor, 
SK. This value is then carried throughout the calculations for this trace. 
If the positional data is included, the first output is a formalized print- 
out of this informational heading. Once this has been executed, the main 
program starts a searching routine to identify various possible features 
before calling the specific sub-routines necessary to analyze this data. 
The program first looks for double-valued temperatures in the vertical 
dimension. If these double-values exist the program calls the inversion 
subroutine, INVERS, to calculate the top and bottom of the inversion and 
its spacial and thermal magnitude. This subroutine first calls a special- 
ized subroutine, ORDERL, that takes all the input data from that trace 
and checks to insure that the data is in correct order from the surface 
down (consecutive readings) and if it is not, reorders the data in the 
correct sequence. Once the data is in correct order, then the program 
does a quick safety check to insure that it does not have a non-inversion, 
(where the temperatures first increases with depth and then decreases, 
giving not an inversion, but a positive surface gradient followed by the 
normal main thermocline, a feature that the main program might identify 
as an inversion). If this condition does exist, the subroutine terminates 


and returns the flow to the main program for further investigations. If 
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an inversion exists ,then, by the use of do-loops, INVERS calculates the 
top, bottom, and magnitudes of the inversion as outlined in Section E.3. 
Once the data has been analyzed and these results are available, the 
subroutine then calls four specialized printing subroutines, PRINT3 
through PRINT6, to print out the data in a useful format. INVERS then 
returns the order flow to the MAIN program for further investigations. 

The MAIN program then investigates the thermal interval, delta-T, 
to insure that if the interval is less than 0.2 degrees a thorough 
investigation for the presence of transients is accomplished. If the 
interval is 0.2 or less, the MAIN routine calls the transient investiga- 
tion routine, TRANS, which proceeds to analyze for the presence of transi- 
ents as outlined in E.1. Upon completing the analysis, the results are 
printed out, if any exist, and the flow returns to the MAIN routine. 

Once returned to the MAIN program, the final investigation is made 
into the possibility of multiple thermoclines. The MAIN program calls 
the multiple thermocline subroutine, MANYT, which utilizes analysis 
procedures outlined in Section E.2. Within this subroutine is the pro- 
vision to go to another subroutine, NONGA, designed to handle the more 
specific case, should only a single main thermocline exist. The NONGA 
subroutine uses the procedure outlined in Section D to find the top, 
center, and bottom of the thermocline. In both subroutines, looking at 
thermoclines, the results are printed out by utilizing the specialized 
printing subroutines, PRINT] and PRINT2. 

A flow chart (simplified) for the entire process is shown in Ap- 
pendix I, and the entire program can be inspected by referring to Ap- 
pendix II. It is worthy of note, that this program will not identify all 
three features, transients, multiple (or single) thermoclines, and in- 


versions, in a single trace as the existance of all three would be very 
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rare. However, the program will look for transients in the upper near- 
isothermal layer and then investigate the possibility of an inversion be- 
low this area, or if the inversion does not exist then look at the thermo- 
cline(s) in the region below the near isothermal layer. 

Also, as an additional feature the program can be called to print 
out a graphic display of the input data (Appendix III), thus giving a 
visual representation of the data that is printed out on the output data 
sheets. This auxiliary feature can be called by inserting the following 


cards in the MAIN program immediately following the card serialized0199. 


CALL SEE13(T,Z,ICOUNT) 
which will call the SEE13 subroutine that utilizes the library pro- 
gram, ''DRAW", to give graphic display of the input data set and fits a 
curve to these points. Before using this call, the user must exercise 
care and carefully read the rules for usage of "DRAW" in subroutine 
"SEE13'' as are outlined in Appendix III. 

Another feature included as part of this thesis package, is a pro- 
visional set of cards (Appendix IV) that will convert the program as 
contained in Appendix II into the slightly different machine language 
required when the program is used on the CDC 6500 model computers, in- 
stead of the IBM 630 models for which this program is written. Before 
using this conversion program two steps must be taken, 

(1) all data and system control cards must be removed 
and the remaining portions of the included program 
become the data deck for 'BCD" 

(2) All output formats that involve printed words or 


combinations of printed words and numbers must be 
converted to "H" format. 
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2 Output from the Computer. 


The output from the computer consists of several possible combina- 
tions of data.. Included are examples of a few of the possible combina- 
tions. All the output from a given trace is printed on a single indivi- 
dual page, headed by the trace number within the data sequence and the 
positional data for that slide (if any). Figure 7 shows one of the sets 
of real data that was taken from some data traces made by Tully in the 
vicinity of Ocean Station "Papa" in 1963. Figure 7 is a hand-drawn en- 
largement of a real test trace containing both transients and a single 
main thermocline, typical of that region, used in testing out this pro- 
gram. Figure 8 is a typical output from the program's analysis of this 
trace and lists the positional data (which was partially available) and 
the date and time of the trace. The program then reveals the existence 
of two transients, by printing out the location of the start and end of 
each transient. It also gives the information on the single major thermo- 
cline. In this case the data was read in from a data dekce after being 
transcribed from the trace onto aperture cards. The data interval was 
taken in this case to be 0.1 degrees. Figure 9 shows the typical print- 
out of a hand-constructed case, also used to test program, centered 
around a thermal inversion (Boston 1966). Figure 10 is the data output 
for this created case. Figure 12 is an example of the printed output 
that would be expected from a trace (see Figure 11) containing multiple 
thermoclines. In this case two major thermoclines exist in the trace. 
For interest, Figure 13 shows the output from the drawing routine called 
from SEE13. It can be seen that computer drawn trace of Figure 13 is an 
adequate enough representation of Figure 12 to enable the user to use 
this type of computer output as a rough aid in insuring correct working 
of the computer program or other rough calculations where a visual presenta- 


tion might be desirable. 
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Fig.7 Ocean Station ‘Papa’ — 


23 April 1963 


a7, 


(after Tully 1963) 
BT Trace 


START OF INVESTIGATIONS IN BATHYTHERMOGRAPH TRACE NUMBER 4 
GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 


LATITUDE e e e e * e ¢ e e > 2 2 e @ e e ® ® 36 ~36 @ 1 NORTH 
LONGITUDE e ° ° e ° e ° e e 6 e o e ° ° e & i 13 -44 e 7 WEST 


DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 


DATE e ® e * e ® e © ® ® e ® ® 6 @ ® e e e 10/23/63 
et ee eke! > o ace 6 « « LOO832 


BATHYTHERMOGRAPH SERIAL NUMBER: 
BOM EAMNO 6.46. 6 66) 6 lees « © 9) ses 0039/7 
NO TEMPERATURE INVERSIONS EXIST IN THIS TRACE, 


TRANSIENT DATA AS COMPUTED BY PROGRAM: 


BEGINNING OF FIRST TRANSIENT. .......-.- Zlae 
ENDVOF TRANSIENT ee. « . ee. + « aMete « 55.0 
BEGINNING OF NEXT TRANSIENT ......2.. .~ 60.0 
END@OP TRANSIENT oe sce ss eg ee sf le 87.2 


MULTIPLE THERMOCLINE DATA AS COMPUTED BY PROGRAM: 


CENTER NUMBER 1 LOCATED AT . ..« « « «© « « 148.0 
CENTER OF THERMOCLINE e e e e e e e e e e e e 149 e 6 
TOP OF THERMOCLINE @ e e e e e e e e e e e e 99 e 3 


BOTTOM OF THERMOCLINE . . . -« « « «© «© © «© «© « 196.8 


END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 
NUMBER OF DATA PAIRS ANALYZED - 18 


Fig. 8 Example of Output for Real Data Curve 
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Fig.9 typical Inversion 
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START OF INVESTIGATIONS IN BATHYTHERMOGRAM TRACE NUMBER 2/7 
GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 


PPI es se lw lt ll ll 00.00 
LONG ITUDE ° ° ° e ® ° e ® @ ® e ° ° ® ® e ® e 00 e 00 


DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 


ee, eg gg she ERM e es ves 00/00/00 
SUE tM, Coe OST. 4 kk e AR ew 000000 


BATHYTHERMOGRAPH SERIAL NUMBER: 
BI OSERTAE NO .%. . 6 se 6 s » CD ew ws 000000 


THE FOLLOWING INDICATES THE COMPUTED TOP OF THE TEMPERATURE 
INVERSION, BY NON-GAUSSIAN THEODS: 


TEMPERATURE INVERSION TOP... ....+.2e-. 66.0 


THE FOLLOWING INDICATES THE COMPUTED BASE OF THE TEMPERATURE 
INVERSION, BY NON-GAUSSIAN METHODS: 


TEMPERATURE INVERSION BASE. .....-e«.-.- 160.1 


THE FOLLOWING FIGURES DEFINE THE MINIMUM TEMPERATURE OF THE 
TEMPERATURE INVERSION: 


MINIMUM TEMPERATURE ° ° e ° ° ° ° e ° ° ® e e 1 1 e O 
DEPTH OF MINIMUM TEMPERATUR Re ae 139.4 


THE FOLLOWING FIGURES DENOTE THE TEMPERATURE AND SPACIAL 
MAGNITUDE OF THE TEMPERATURE INVERSION: 


THERMAL MAGNITUDE ° ° e e e ° ° ° ° ° ° ° ° 4 e 0 
SPACIAL MAGNITUDE a ee ee 94.1 


NO THERMAL TRANSIENTS EXIST IN THIS TRACE, 
NO MULTIPLE THERMOCLINES EXIST IN THIS TRACE. 
END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 


NUMBER OF DATA PAIRS ANALYZED - 


Fig. 10 Example of Output for Inversion 
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Fig.l! Hand drawn Bathythemogram 
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START OF INVESTIGATIONS IN BATHYTHERMOGRAM TRACE NUMBER 
GEOGRAPHIC POSITION OF BATHYTHERMOGRAPH TRACE: 


LAT ITUDE e e e e e e e e e ® e ® e e e e e e 
LONG ITUDE * e € 9 ® ® @ e ° e e @ @ © 2 e @ ® 


DATE AND TIME OF BATHYTHERMOGRAPH TRACE: 


DATE e e 9 e o e e e ® ° e e e e e ® @ e e e 
TIME @ e e e e e * e e e e e e e e e e e e @ 


BATHYTHERMOGRAPH SERIAL NUMBER: 
Bi wSermaeNO. . 5 6 4% wes ws © 6 Soe ie ees 

NO TEMPERATURE INVERSION EXIST IN THIS TRACE, 

NO THERMAL TRANSIENTS EXIST IN THIS TRACE. 

MULTIPLE THERMOVLINE DATA AS COMPUTED BY PROGRAM: 
CENTER NUMBER 1 LOCATED AT ........ @ 
CENTER NUMBER 2 LOCATED AT .... 2. e2e- 
CENTER OF GAUSSIAN THERMOCLINE ....... 
TOP OF GAUSSIAN THERMOCLINE .......24.. 
BOTTOM OF GAUSSIAN THERMOCLINE ....... 
SEPERATION #1 BETWEEN THERMOCLINES ..... 
CENTER OF GAUSSIAN THERMOCLINE ....... 
TOP OF GAUSSIAN THERMOCLINE ........ . 

BOTTOM OF GAUSSIAN THERMOCLINE .....5 . ee. 


END OF DATA FOR THIS BATHYTHERMOGRAPH SLIDE. 
NUMBER OF DATA PAIRS ANALYZED - 


ke 


000.00 NORTH 
000.00 WEST 


00/00/00 
00000 


000000 


72.0 
144.0 
75.0 
Dom) 
94.4 
100.0 
138.0 
104.9 


15250 


Fig. 12 Example of Output for Multiple Thermoclines 
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All the preceeding traces were read in in data deck format, with 
varied intervals of delta-T, 0.1 degrees in the first example, 1.0 
degrees in the second case, and 0.5 degrees in the latter case. There 
are other possible combinations, but these three cases outline the general 
format and the answers that can be expected. The other cases are just 


various combinations or omissions of the above three cases. 
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6. Discussion and Conclusions Concerning Results. 


Although no error analysis was undertaken for this computer pro- 
cess, it is obvious that there are some errors built into the program 
that could perhaps be reduced by further investigation and refinement. 
In several segments we use a Gaussian based method and thus there is 
some error introduced in the computation of the top and bottom of the 
thermocline due to having to compute a sigma in the program. Also, in 
this program the multiplication factor, SK, is calculated only roughly 
based on the thermal interval. This factor can be improved, but cannot 
be exact in the computer, due to truncation error; besides, much storage 
space is needed to store the large table of values applicable where each 
tabular figure is only good for one specific set of data, 

There are certain errors inherent in the computer, but since we do 
not require great precision, these errors can be ignored due to the fact 
that the bathythermograph itself is not a precise instrument. These 
errors include truncation and round-off errors in every computation. 
Since we are requiring printout to saat one decimal place it is doubt- 
ful that these errors will ever be noticable. Also there are errors 
built into the functions that are interral to the computer, such as the 
"square root," "multiplication," and "addition" routines, but again these 
will not be significant in the accuracy that we require. 

There are certain advantages to the use of data tapes versus magne- 
tic cards, when great volumes of data are to be processed. The time 
required to punch the data cards and to check the results is a function 
of the experience of the operator. The computer time is approximately 
25 seconds for card input computations, plus additional time if the draw- 


ing of graphs is required (approximately three minutes per graph). This 
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time can be cut approximately in half if the program compiled by the 
computer is put on tape and data input is in that form, as the compiler 
takes approximately half the time now used in making a computer run. 

There are other equipment errors, that as yet have not been evalu- 
ated,that are concerned with the input data, and are still vital to the 
person desiring to analyze the traces by use of this program. Such items 
as hysteresis, where a double trace is formed on the bathythermograph 
slide, due to slight differences in the temperature when raising the 
Bathythermograph as compared to when lowering it cause inaccuracies , 
dull styluses can give traces that are too wide to be useful. Damaged 
equipment that gives unacceptable temperature fluctuations, such as 
variations from bucket thermometer temperatures, or other signs that the 
bathythermogram is inaccurate although the trace itself appears accept- 
able is another input error. Also there are errors inherent in the digiti- 
zer that can vary depending upon the equipment used and the usual problems 
inherent to electronic curve followers or other devices. 

The results of this program, in the fifty test cases run, appear to 
give slight errors as compared to the hand analysis of the same data. 

Ten of the test cases were hand constructed curves, designed to test 

the program for specific features or combinations of these features, and 
the remaining data were real data obtained from the Fleet Numerical Weather 
Central or taken from other written works, such as Tully and Dodimead 
(1957) and Pattullo and Cochrane (1951), and a few traces taken in 

Monterey Bay. Both tape inputs from Fleet Numerical Weather Central 

and punched tape data decks were used. There were no apparent differ- 
ences between the outputs from the digitized data and the data on 

aperature cards, as might be expected assuming that care was taken by 


the person picking off data points from the curves. The largest error 
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in the tested cases, was found in a very steep thermocline, real data 
case, where the depth error was 4.5 feet for the top of the thermocline 
and 3.0 feet for the base of the thermocline in a thermocline that was 76 
feet in vertical magnitude. In this case the center of the thermocline 
was only 0.3 feet from the hand calculated value. The thermal interval 
was 1.0 degree Farenheit and the thermocline in question was analyzed 

by the Gaussian technique. 

Thus it can be said that this program accomplishes, with acceptable 
accuracy, considering the analytic methods used, in seconds what a well 
trained operator would take tens of minutes to accomplish. 

This program can be expanded, by the inclusion of additional sub- 
routines, and their respective 'call' cards in the MAIN routine without 


greatly increasing the time needed in computing the desired results. 
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ae Recommendations. 

As previously stated, there are other features that can be added 
to this program by merely including additional subroutines. Additionally, 
this same programming technique can be used, with a few basic changes, to 
analyze other oceanographic data curves. Before accomplishing any con- 
versions of this program to handle other data, the user must first look 
at the basic form of the traces he plans to analyze. Traces, other than 
temperature, may not be describable using the same parameters as des- 
cribed by Tully (1964). Such items as the magnitude of the fluctuations 
of other quantities must be examined as they may exceed the number that 
this program can handle. For example, it may be common, in some other 
oceanographic parameter, to have more than three stable features corres- 
ponding to the thermoclines in this program, and thus the program would 
have to be modified to allow for additional curve computations. Another 
correction that would have to be made would be the rewording of the print- 
out statements. 

Further investigations should be undertaken, to include the improve- 
ment in accuracy of sigma and SK parameters and the possibility that such 
curves as the hyperbolic tangent or Chi-squared curves might give a better 
Curve tit. 

As this program has yet to be run on the CDC 6500 computer of the 
Fleet Numerical Weather Central, this should be undertaken as soon as 
possible. For this type of program, tne speed of computation on the CDC 
6500 computer is approximately three times that of the IBM 360. This 
CDC 6500 computer also contains more storage space than the IBM 360 
computer and thus longer tapes or card inputs may be used. An additional 


advantage of this increased storage is that investigation could be shifted 
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from the present two dimensional set up (temperature versus depth) to 
a three dimensional investigation (multiple traces analyzed simultaneously) 


giving a field representation instead of just specific points. 
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Initialize NBT 


PRINT OUT 
Print out "Start of investigations into 


trace characteristics for bathythermo- 
graph trace number" 





CALL SEE13 READ: T(T) ,2(1) 


Does T(I) equal zero and simultaneously 
Z(1I) greater than: zero? 


YES 
Does T(I) equal zero and simultaneously 
Z(1) equal zero? 











Increment counter 








READ: 


LAT,XLAT, LONG, YLONG, QUADI1, OUAD2, 
DATE1, DATE2, DATE3, TIME, SERIAL, to 
determine the time, serial number, and 
location of bathythermograph drop 








4) Call PRINT? (T) 


Test for repeat values of temperature 
in individual trace. 


a Call INVERS (A) 





uw 
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PRINT OUT 


Print out, "'No thermal transients exist 
in this trace." 


PRINT OUT 


exist in this trace." 













Test thermal interval to see whether or not 
it is worthwhile investigating for the 
resence of transients. 


(2) CALL TRANS CK) 


Soe 






Investigate the data set to see whether there 
exist more than one major thermocline. If 
more than one thermocline exists, each 
thermocline is investigated as to the location 
of its center, top, and base. If only one 
thermocline exists then a more complete analys 
Sis ig made of this single thermocline. 


£3} Call MANYT (me) 


PRINT OUT 
Print out, "Completion of all investiga- 
tions into bathythermograph trace number." 







Print out, "No multiple thermoclines 
} 






NO 








WV 
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Is ICOUNT greater O ° 


WY 
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Subroutine PRINT7 | 


Inputs from: MAIN 
LAT, XLAT, LONG, YLONG, QUAD1, OUAD2, 
DATE1, DATE2, DATE3, TIME, SERIAL 


[0 


QUAD2=1 NO 


YES 


QUAD1=1 





QUAD2=1 








PRINT OUT 


Print out; 'NORTH®™ latitude and "EAST" 
y | longitude and also date, time, and serial 
of bathythermograph. 


PRINT OUT 


Print out; "NORTH" latitude and "WEST" 
longitude and also date, time, and serial 
of bathythermograph. 


PRINT OUT 


Print out; "SOUTH" latitude and "EAST" 
longitude and also date, time, and serial 
of bathythermograph. 


oY 









PRINT OUT 


Print out; "SOUTH" latitude and "WEST" 
longitude and also date, time, and serial 
of bathythermograph. 
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Subroutine INVERS 


Inputs from: MAIN 
Z(I), TCI), ICOUNT 


Define SK in terms of temperature interval. 
Define and set do-loop counters. 


ae fe 
Call ORDER2 © 


Test entire data set to insure inversion 
exists by checking sign changes of: 
T(NA1) - (TC(N)). 


By do-loop, compute first difference, A(N), 
on Temperature, and first difference, B(N), 
on depth, terminating loop if sign of 

A(N) changes 


By do-loop, compute second difference, 
G(L), based on B(N). 







Set up initial testing loop and set 
TEST1L equal to G(L). Run loop until 
sign of TEST1 changes and note loop 






Compute top of inversion, ZCl, based on 
TEST1L termination counter, and also comp- 
ute temperature at top of inversion, 

EME. 







ao 


SS a TS ee Se eg ES A rE a ES. plese 


© 


By do-loop, compute first difference, 
C(131), on temperature, and first differ- 


ence, D(I31), on depth, terminating loop 
pn of C(I31) changes. 


Call ORDER} 


By do-loop, compute second difference, 
P(I132), based on D(I31). 


Compute spacial magnitude of inversion, 
AUPP, based on loop counter NNCO, 





By do-loop, compute first difference, 
E(133), on temperature, and first differ- 
ence, F(133), on depth, terminating loop 
at counter, based on length of data set, 
ICOUNT. 


By do-loop, compute second difference, 
R(1I34) on F(I33). 


Set up second testing loop and set 
TEST2 equal to R(I35). Run loop until 
sign of TEST2 changes and note loop 
counter at this point. 


Compute base of inversion, ZC2, based 
on counter when TEST2 sign change 


encountered. 
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eo 


Compute minimum temperature, TMIN, and 
its corresponding depth, ZMIN, based on 
ounter in R(134) do-loop. 


eit as _— F ) 
[an mm  [-—*6) 
CH) 
[esti renne sf 


RETURN 
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Subroutine ORDER2 


is 


Inputs from: INVERS 
X(I), III, M 


This subroutine takes the input data and 
rearranges it, if necessary, in correct 
descending order, by using a double do- 
loop. 


DO I =M, III-3 


IP2 = lf 1 


DO I = IPl, III-3 


YES 


SUBST = X(I) 


X(I) = X(J) 


X(J) = SUBST 


CONTINUE 


rare 


RETURN 


OZ 





| Subroutine ORDER] | 
[=eee 


= Se ES ERE > aE eS 









Inputs from: INVERS 
X(1), T(1), III 


mle 
oo Se, =e a eyo - a ss oe oe Soi RES) ee 


This subroutine, similar to ORDER2, takes the 
input depths and temperatures and sets them in 
correct order using a double do-loop. 








ae baat | 
SUBST = X(I) | 


oe 
——_ ~ a 
| 


| X(I) = X(J) J 


Ss 





oe Oe Ce ee 6 + re ome ae es me re re ee ne ee = 


| SUBST = T(I) 


= a 
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TCI) = T(J) 


CONT INUE 
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Subroutine PRINT3 


Inputs from: INVERS 
ZCl 


PRINT OUT 
Print out, ‘top of temperature inversion 


based on on-gaussian methods" 





RETURN 
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Subroutine PRINT4 


ic 


Inputs from: INVERS 
“CZ 


PRINT OUT 
Print out, "Base of temperature inversion 


based on non-gaussian methods." 





RETURN 
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Subroutine PRINT5 _J 






Inputs from: ravers 
TMIN, ZMIN 





PRINT OUT 
Print out, ‘Minimum temperature in tempera- 
ture inversion" and 'Minimum temperature 
found at depth." 


po 





——— Se 


RETURN 
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Subroutine PRINTO6 


Inputs from: INVERS 





ZUPP, TUPP 









Print out, "Spacial magnitude of temperature 
inversion" and "Thermal magnitude of tempera- 
ture inversion." 
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Subroutine TRANS 


Inputs from: MAIN 
Z(1), TCI), ICOUNT 


Set do-loop counters, based on size of 
data set, ICOUNT. 


Bydo-loop compute first difference, 
Z6(N), on depth. 









By do-loop, compute second difference, 
2/7 €14)) on Zon 






By loop, check sign of second difference 


a to see if it is greater than zero. 





ee ee ee 






PRINT OUT 


Print out, ‘Beginning of first transient" 





——— aw ee ee 


YES NN.LE.M 


CONTINUE 


YES 


NN.LE.M 
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Increment NN | 


Z7(NN).gt.0? 










| YES 


YES = = ' 
3 


i NN 


INCREMENT NN NN ES 








YES 





PRINT OUT | 


Print out, "End of transient." | 





PRINT OUT 


Print out, "Beginning of next transient" 







NN. LE.M? 
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Subroutine MANYT 


Inputs from: MAIN 
Z(I), T(I), ICOUNT 











Set do-loop counters based on size of 
data set as inputed. 


By do-loop, compute first difference, 
Z6(N2), on depth. 


By do-loop, compute second difference, 
Z8(N1), on Z6(N2). 


Initialize counters 
By do-loop, check all values of Z8(N1) 
to see if greater than zero. 


Increment counter 


PRINT OUT 


Print out "Center number ... located at" 





Set value to 
BIGl 


By do-loop, compute the largest value of 
Z6(n) in thermocline. Set BIGDP1 equal to 
0.5 (Z6(max) + Z6(max + 1)). 







Call GAUSS3 i) 


PRINT OUT 
Print out, "seperation #1 between thermo- 


clines." 





Reinitialize counters 
Set value to: 
BIG2 









Eye o een Coma the largest value of 
Z6(N) in the next thermocline. Set 


BIGDP1 a to 0.5(Z6(max) + Z6(max + 1)). 


cait causss fe ——*{ R) 






PRINT OUT 
Print out, ‘Seperation #2 between thermo- 
clines." 
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YE 






Reinitialize counters | 
CALL GAUSS3 i —"( Rr) 


7 i aca een of thermoclines fen greater than 
If there is only one thermocline in the data ae 
set, then call NONGA. 


| PRINT OUT 
Print out, "End of data for this bathythermo- 
gram trace'’ and "Number of data pairs 
analyzed." 















RETURN 


(s 
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Subroutine GAUSS3 


“ 


Inputs from: MANYT 
Z(I), TCI), ICOUNT, NQN1 


Compute value of SK, based on the tempera- 
ture interval. 


Call ORDER3 ae 


Compute Z1(N) = (Z(N) + Z(n+1)/2.0) and 
Z2(N) = Z1(N) x Z1(N). 


Sum: Z1(N) and Z2(N) within do-loop. 


Compute gaussian center of thermocline, 
ZC = (1.0/(total point - 1))xsum Z1(N). 


PRINT OUT 
Print out "Center of Gaussian thermocline 


is located at depth." 





Compute: ZAVG = (1.0/(total -1) x sum 
of Z2(N); SIGMA = ZAVG - ZCxZC; 


Compute top of Gaussian thermocline 
= = x_ SIGMA 


Wy 
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PRINT OUT 






Print out "Top of gaussian thermocline is 
located at depth" 


Compute base of gaussian thermocline, 
ZBT = ZC + SK X SIGMA 


PRINT OUT 


Print out "Base of gaussian thermocline is 
located at depth" 





Compute magnitude of thermocline, from 
WW = ZBT - ZT. 


PRINT OUT 


Print out "Magnitude of thermocline." 





RETURN 
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© 


Subroutine ORDER3 


Inputs from: NONGA 
ZNT ) 5 KT 


This subroutine arranges the input data, 
if necessary into correct descending order 
from the surface down. 


a a ae re em eee © ee © me en 9 ee Oe a 


DO I = 1, KT-1 
= IPl =1+1 
DO J = IPl, kT-1 









YES 


X(J) = SUBST 


CONTINUE 


Subroutine NONGA 


Inputs from: MANYT 
ZCI )e 1 (1). ICOUNT 






Compute SK, based on the temperature 
interval. 


Compute first difference, by do-loop, 
Z1(N), on depth. 


Compute by do-loop, second difference, 


Z2(n) sone Zl CNe 












Set up testing loop, TEST(KK) equal to 
Z2(KK) and execute loop until there is 
a change in sign of TEST (KK). 






Compute the center of the thermocline, 
ZC, based on the counter number when 
TEST loop terminated. 


Set up localization loop to find extent of 
the upper portion of the thermocline 


Call GAUSS1 (AA 


Set up localization loop to find extent of 
the lower portion of the thermocline 
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Call GAUSS2 
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Subroutine GAUSS1 | 


2 a aaa RS pe Coan 


Inputs from: NONGA 
Z2(I), NT,  ICOUNT 
Compute value of SK, based on input 
value of ICOUNT. 
| Compute; Z1(N) = (Z(N) + Z(N41))/2.0 | 


| CBee: Z2(N) = Z1(n) x Z1(N) | 


Increment N 



















































eel 


By do-loop, compute SUM1 and SUM2, based 
on sums of Z1(N) and 2Z2(N). 


| Compute; 2C = (1.0/total -1) x (SUM1), the | 


center of the gaussian thermocline. 


a 


| Compute; ZAVG = (1.0/total - 1)x(SUM2); | 














SIGSQR = ZAVG - (ZCxZ2C); ZT = ZC + SIGSQR 







Call PRINT1 


RETURN 
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Subroutine GAUSS2 


Inputs from: NONGA 
Z(1),>_ Ni, ICOUNR 


Compute value of SK, based on input value 
of ICOUNT. 


wees eee ee 














owe ae ee by 


| Computemzicn) = (zea + Z(N+1))/2.0 


| Compute Z2(N) - (Z1(N) x Z1(N)) | 


Increment N 


By do-loop, compute SUM1 and SUM2 based 
on sumes of Z21(N) and Z2(N). 
| Compute; ZC= (1.0/total - 1) x SUM1 


Compute; ZAVG = (1.0/total -1) x SUM2; 
SIGSQR = ZAVG + (ZC x AC 


Compute base of thermocline, ZBT = ZC + 
SK x SIGSQR. a a 
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ES AAD ED Ce 


Subroutine PRINT1 | 


-- a SS NS a RL Ee 


Inputs from; GAUSS] 
ZT, 826, ICOUNT 


PRINT OUT 








"Number of data pairs utilized" 


PRINT OUT 


Print out "Center of thermocline is 
located at depth. 


PRINT OUT 


Print out ''Top of thermocline is 
located at depth." 


RETURN 
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oe See ie OSE 


{ 


ae 


Subroutine PRINT2 4 
(2 


a a 


} 


———-1L---—., 


; Inputs from: GAUSS? | 


ie a aac 


ee me ee ee eee ee eee 


PRINT OUT 


Print out ''Base of thermocline located at 
depth." 


me eee em eee ee ee 


_4 
or) 


_ 
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INPUTS TO THE MAIN PROGRAM: 


ARE THE BT TRACE DATA INPUTS 


THE FOLLOWING 


ae READINGS IN DEGREES 


A 
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EMPER 
ARENH 


b— L. 


TET) 


ORRESP 
OR EGE 


CORRESPONDING TD TETD&S 
EET. 
ARE THE POSITIONAL DATA (SEE USER'S NOTF-3 


Z(1) 


OLLOWING 


THE F 
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WHOLE DEGREES OF LATITUDE 
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AND HUNDREOTHS OF MINUTES OF DEGREE 


XLAT 


WHOLE DEGREES OF LONGITUDE 


LONG 
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ORNER. 


SURFACE DOWNWARD, 


INSURES CORRECT, 


ORDER2 


OUVUOVU OU UN OOU0 UO 


PRINTS OUT OEPTH OF TOP OF INVERSION 


PRINT3 
PRINTS 
PRINTS 


PRINTS OUT DEPTH OF BASE OF INVERSION 
- PRINTS OUT MINIMUM TEMPERATURE AND ITS DEPTH. 


PRINT6 - Ae OUT SPACIAL AND THERMAL MAGNITUDE OF THE 
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